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Abstract 
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I. INTRODUCTION 



The accelerating expansion of the universe has been estabhshed by several cosmological 
observations, such as those from type la supernovae Hi, cosmic microwave background 
radiation js, 4|, large scale structure ^| and weak lensing [gI. There are two ways to explain 
this phenomenon. One is to add dark energy to modify matter and the other one is to modify 
gravity in the Einstein's equation. The simplest version for the latter is f{R) gravity 7H9|, 
which is by extending the Ricci scalar of i? to a function of f{R) in the Einstein-Hilbert 
action. As a result, the late time accelerating universe can be realized in f{R) gravity. Many 
viable f{R) gravity models have been constructed by satisfying various conditions as well 



as constraints from cosmological observations 



9| 



The finite-time singularity problems 10| have been examined in many modified theories 



such as f{R) models 
Horava-Lifshi tz g ravity 



modified Gauss-Bonnet models 3, f{T) models [l3|, modified 



out in Refs. 



16 



14j and non-local gravity models 15|]. Recently, it has been pointed 



17l | that some of the viable f{R) models contain one kind of the finite- 



time singularities, leading to a divergence of curvature, but it can be avoided by taking a 
fine-tune initial condition 18, lol. However, this kind of singularities must be induced and 
can not be avoided when the local background density of matter becomes dense 20| . Even 
though the singularity depends on the background density as well as the model parameters, 
it happens in a finite time. This behavior could exist in many physical systems, such as 
cluster, galaxy, nebula collide, and star collapse. However, if an additional i?" term with 



1 < n < 2 is introduced in the viable f{R) models 2ll-l23|. the singularity can be avoided. 



We note that adding the R^ term to the viable f{R) gravity models could also lead to the 
unification of dark energy with inflation j^]. We also remark that the curvature singularity 
arises naturally in the viable f{R) models unless some fine-tuning is taken 23|. 

In this paper, we first show the singularity problem in the popular viable f{R) models and 
then try to modify them to remove the singularity. We will explore the possible cosmological 
consequences in these modified models under the Minkowski background. In particular. 



24 



251, which is a characteristic 



we will analyze the scalar mode of gravitational waves 
signature to distinguish f{R) gravity from general relativity (GR). We will also show that 
the Minkowski approach still holds when the curvature close to the singularity. We use 
the natural unit c = h = 1 with Mpi = G~^/^ ~ 1.2 x lO^^GeV and the metric gfj_^ = 



diag{-, +, +, +). 

The paper is organized as follows. In Sec. [Tll we study the curvature oscillation in the 
viable f{R) gravity models with a highly dense background. In Sec. IIIIl we include the 
term in the models to prevent the singularity problem. The scalar mode of gravitational 
waves is also discussed. The conclusions are given in Sec. IIVI 



II. CURVATURE OSCILLATION IN f{R) GRAVITY WITH HIGHLY DENSE 
BACKGROUND 

The action of f{R) gravity is given by 

5 = ^1 d^x^{R + f{R)) + SU9,u, V) , (2.1) 

where g is the determinant of the metric tensor g^^, Sm is the action of matter, = SttG = 
M~i^, T denotes the matter field with the minimal coupling to gravity, and f{R) is an 
arbitrary function of the Ricci scalar R. By varying the action (12. ip with respect to g^j,^, we 
obtain equation of motion 

(1 + fR)R^.u -liR + f)9,.v + {g^^uU - V^V.)/ = K^T^,, (2.2) 

where / = /(-R), the subscript R denotes the derivative with respect to i?, i.e., = df /dR, 
is the covariant derivative, □ = g^''^V is the d'Alembertian operator, and T^i, is the 
matter energy-momentum tensor from the local matter distribution. 



A. Viable f{R) Models 

It has been widely accepted that a viable f{R) gravity model has to satisfy the following 
conditions js]: (a) 1 + //j > for i? > i?o, which keeps the positivity of the effective 
gravitational coupling and avoids anti-gravity, where Rq is the present background curvature; 

(b) fnR > for R > Rq, which gives the stability condition of cosmological perturbations; 

(c) f{R) — )■ -R — 2A in the large curvature regime {R ^ Rq), which realizes the ACDM 
behavior at i? ^ Rq; (d) a stable late-time de-Sitter point; and (e) passing local gravity 



constraints, including those from the equivalence princip 
conditions, many viable f{R) models have been proposed 



e and solar system. Under these 
9| . In Table HI we give the explicit 
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forms of the popular viable /CR) models in the literature |9|, where (i), (ii), (iii), (iv) and 
v) correspond to Hu-Sawicki 26], Starobinsky Tsujikawa 28], the exponential gravity 
29l-l33| and Appleby-Battye (AB) [3J, l35| models, respectively. 



TABLE I. Explicit forms of f{R) in (i) Hu-Sawicki, (ii) Starobinsky, (iii) Tsujikawa, (iv) the 
exponential gravity, and (v) AB viable models 



model 


f{R) 


Constant parameters 


(i) 


ciWRhs)" 
^HSe2(i?/j?Hs)''+l 


ci, C2, p{> 0), Rus{> 0) 


(ii) 


-XRs l-(l + |J)"" 


A(> 0), n(> 0), Rs 


(iii) 


-fiRT tanh(^;^) 


fi{> 0), Rt{> 0) 


(iv) 


-I3Re (1 - e-^/^E) 


P, Re 


(v) 


-gR + gRABln\-'^^^^r-'^' 


9, b, Rab 



B. Curvature Oscillation in f{R) Gravity 

We start from the trace of the field equation in Eq. (I2.2p . given by 

RfR-2f-R + 3DfR = K^T, (2.3) 

where T = g'^'^T^^ is the trace of the energy-momentum tensor. This equation is reduced to 
GR {R = — K^T) if f{R) = 0. Clearly, it contains an extra degree of freedom beyond GR 
when //J 7^ 0. Note that Eq. fl2.3p is a fourth order field equation in comparison with the 
second order one in GR. This fourth order equation also leads to the oscillation behavior of 
the Ricci scalar. 

Before calculating the curvature oscillation in the viable /(-R) models, we briefly introduce 
the framework of our study. It is known that the curvature in the dense matter regime, such 
as inner-galaxy, nebula and star collapse, is much bigger than the background curvature of 
the universe. Note that the average density of the universe and the density inside the inner 
galaxy are 10~^^ and 10~'^^g/cm^, respectively. In the large curvature limit, the viable f{R) 
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models in Tabl^ can be reduced into power law and exponential types, given by 



-XR, 



ch 



R. 



ch 



R 



2n 



ill): /(/2)^-A/2,,(l-e-^/«-) 



(2.4) 
(2.5) 



referred to as Type-I and II, respectively, where Rch is a dimension-two constant. The Hu- 
Sawicki and Starobinsky models belong to Type-I in Eq. (12. 4p . while Tsujikawa, AB and 
exponential models correspond to Type-II in Eq. (12.51) . In order to explain the accelerating 
universe, the parameter R^h should be the same order as the background value of Rq at the 
current epoch. 



1. Curvature Oscillation in Type-I f{R) Gravity 

We consider f{R) gravity in dense, locally homogeneous and isotropic perfect fluid with 
the non-relativistic matter density p^, and assume that the density changes homogeneously 
in time and is much denser than the critical density pc- Then, the trace of the energy- 
momentum tensor can be expressed as 

T = -To + (2.6) 



tch 

where Tq = pm — 3Pm •* ~ pm is the initial background density and tch is the characteristic 
time. By defining a dimensionless variable 

u = Rch/R, (2.7) 

/ and ffi of Type-I in Eq. (l2.4p become 



/ - -XR 



ch 



1 - 


( Rch 


V R 


Rch^ 


2n+l 


R J 


) 



2n 



-\Rch (1 - n^") , (2.8) 



/k~ -2n\ \-^\ = -2nAM^"+\ (2.9) 



respectively. The trace of the field equation in Eq. (12. 3p results in the oscillation behavior 
in the large density limit, given by 

Dfn = 2nA— (u^'^+i) = 2n {2n + 1) A {u^'^u + 2nu^''~^u'^) 

= -I j/t'To (l + - Rch [2n\u^^ - 2A (1 - n^") + , | (2.10) 



where we have assumed that the curvature and energy density depend on time only. By 



using the dense background h?Tq/R, 
simphfied to 



u + 2n h 



u 



ch - R/Rch{= u 

K io 



Pm/pc > 1, Eq. fl2A0|) can be 



u 2n{2n + l)X 



t 



tch 



ch 



3u 



(2.11) 



With the rescaled variables y = (5u and r = 7 ^t, we obtain 



y" + 2n 



12 



-2n 



1 + 



r \ 



- 



Teh J 



(2.12) 



where 
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Rch 

2 Qn {2n + 1) A 



R. 



ch 



R, 



ch 



2n+2 



(2.13) 
(2.14) 



and the prime denotes the derivative with respect to r. Note that /3 is a dimensionless 
parameter, which rescales y = (3u = k'^Tq/R to be unity when the background density 
is stationary at the initial value of i? ^ /^"^Tq, while r is also a dimensionless variable, 
related to the physical time by a constant and time dimension factor 7. Since tch = I'Tch, 
7 can be estimated under various backgrounds and model parameters. For example, 7 ~ 
(400,4 X 10-3)s with n = (2,3) and p„ ~ 10~'^^g/cm^. 

In Fig. [H we show the evolution oi y = k^Tq/R as functions of the rescaled time r in 
the Type-I models with initial conditions of = 1 ^^'i^ Vo = 0- From the figures, we see 
that the curvature singularity of i? — ?■ 00 appears when y = h?Tq/R reaches zero. Clearly, 
the singularity exists in a time shorter than the age of the universe. The evolution equation 
(I2.12P only depends on the characteristic time Tch-, so that the divergence always happens at 
the order of Tch- 



2. Curvature Oscillation in Type-II f(R) Gravity 

We now consider the divergent behavior in the Type-II viable f{R) models. In the dense 
background density, from Eq. (12.51) we have 

/^^-Ae-«/^-, (2.15) 




FIG. 1. Evolution of y = k?Tq/R as functions of the rescaled time r in the Type-I models with 
n=2 and (a) = 50 and (b) Tch = 100, where the dashed lines correspond to R = —k?T = 

K2ro(l+ r/Tefe). 



Substituting Eqs. (12. 5p and fl2.15p into the trace equation (12. 3p . we obtain 



X — X + 



xR, 



ch 



3A 



'"'^"-e^ ( 1 + _ 1 ~ 0, 



3A 



t 

tch 



(2.16) 



where x 



u 



R/Rch- Because the e-folding contains a variable x, it is hard to rescale 



this equation into a background independent equation as that in Eq. (I2.12p . However, we 
can stih redefine some parameters to modify the evolution equation to 



/2 



y" + i(3y-'-2f- + y'ef^/y 

y 



-) 

Tch J 



(2.17) 



where /3 = ti^TQ/Rch, y = f3x~^ = li^To/R and the prime denotes the derivative with respect 
to r, defined by r = C,~^t with = 3X/Rch- The initial conditions can be determined easily 
to be ?/o = 1 ^-iid y'o = ^ when the background is stationary. The time scaling factor can 
be estimated as ^ ~ 4.3 x lO^^s. In Fig. [21 we illustrate the evolution of y as functions of r 
in the Type-II viable f{R) models. The figures show that the divergent behavior depends 
very strongly on the background, and the singularity appears not only in a finite, but a very 
short time. 
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FIG. 2. Evolution of y = k'^Tq/R as functions of the rescaled time r in the Type-II models with 
'Tch = 0.1 and (a) 13 = 16 and (b) (3 = 18, where the dashed lines correspond to ii = —k^T = 

K^Toil+T/Tch). 



III. PREVENTING THE SINGULARITY PROBLEM 



A. The Term in f{R) Gravity 



It has been shown that the singularity could be prevented by an additional / M'^^'^~^'> 
term with 1 < m < 2 j23)|. Since the inflationary 36|-|39[ and evolutionary 40[ model of 
f{R) = R + R'^/M'^ has been well-discussed, it is reasonable to examine its behavior with 
m = 2. In the large curvature regime, we rewrite the viable f{R) models in Eqs. f l2.4p and 
(12. 5p plus the R^ term as 



(/) : f{R) = fnss + fn^ 



-XR, 



ch 



R. 



ch 



R 



{II) : f{R) = fE.p + fR2 ~ -XRck (1 - e~^/^- 



2n 



) + 



R^ 



R^ 
M2 ' 



(3.1) 



respectively. 

We now examine whether the singularity problem can be resolved in the above modifled 
viable f{R) models. In the modifled Type-I models, with the similar procedure in Eq. (12.111) . 
the trace equation (12. 3p can be rewritten as 



12 



2n h giy 

y 



-(2n+2) I ylf 



2y 



12 



-2n 



Teh J 



0, 



(3.2) 





FIG. 3. Legend is the same as Fig. [T]but with r^h = 50 and (a) gj = I and (b) gj = 10 



where y = k^Tq/R, t = and gj arises from the addition term, given by 

2n+2 

It J I K 1q \ 

9i 



(3.3) 



Xn{2n + 1)M2 V Ri 

with the time rescahng factor 7^ = (6An(2n + 1)/Ri){Ri/{k'^To)Y''^'^ and Rj = Rch- In 
Eq. fl3.2p . the curvature singularity could be prevented by the additional factor gj since 
there is a huge restoration force in a large curvature regime. Similar to that in Eq. (12.121) . 
one finds that i? — )■ 00 if y — in Eq. (13. 2p . Moreover, Eq. (13. 2p is also a scale independent 
evolution equation with the period of the oscillation depending on the time scaling factor 7. 
We illustrate the oscillation behavior in Fig. [31 Clearly, the singularity can be avoided when 
we include the R^ term in the Lagrangian. It still holds even when ^ 1, but the amplitude 
is strongly related to the dimensionless positive parameter gj in Eq. (13. 2p . We note that the 
singularity appears at r ^ 40 if there is no R^ term. We also note that Rj is determined 
from cosmological constraints, which should be the same order as the cosmo logical constant. 



In the Type-II f{R) models, the evolution equation is 
./2 / 2y'^ 



y" + [Py-' - 2) ^ 
y 



Que' 



y 



y 



1 + 

Teh) 



0, (3.4) 



where y = k'^Tq/R, t = ^ ^t, ^"^ = 3A/i?// and gjj = ^^e^- Similarly, as in the modified 



Type-I models, the singularity could be eliminated when the R^ term is added. In Fig. HJ 
we show the curvature oscillation with two different values of gu. Clearly, its amplitude of 
the oscillation still strongly depends on gjj. 
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'0123456 0123456 
T xlO"* X xlO"* 

FIG. 4. Legend is the same as Fig. [2] but for = 10~^ and (a) gjj = 1 with /? = 25 and (b) 
gjj = 10-2 with /3 ~ 20.4 . 

From the above results, one can easily conclude that the oscillation behavior is deter- 
mined by the dimensionless constants Qui- Consequently, if the curvature is large enough, 
the curvature oscillates as a simple harmonic oscillator with a driving term "k^T". If the 
curvature is not large enough, it is dominated by the original parts of the f{R) models, 
whereas the singularity is removed by adding the term in the large curvature regime. 



B. Cosmological Phenomena in Modified Viable Models with R"^ 

Although we can avoid the singularity with introducing the term in the Lagrangian, 
the oscillation behavior still exists. This behavior might appear in some physical systems 
such as the scalar mode of gravitational waves, which has been recently discussed in Ref. {25 1 
for the viable f{R) models. The graviton in GR is a spin- two massless particle with two 
spin polarizations, corresponding "plus" and "cross" modes, respectively. The scalar mode 
is an extra mode of gravitational waves, coming from the additional degree of freedom 
in f{R) gravity and the non- vanishing trace equation in vacuum (we can estimate it by 
using Eq. (12. 3 p with T = 0). The scalar mode of gravitational waves propagates in the 
vacuum like a massive scalar field nhf = m^hf, where Rmin is the background curvature, 
^21 (l+MR^ _ ^ \ and hf = —j% — . = ^Tr/Ti^'"r^ ■ Because of the same 
origin, the density increasing system could be a source of the scalar mode of gravitational 
waves when the the curvature oscillation amplitude {6R/ Rmin) is large. 
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3^/ "Rm in ( a ) 'R/ "Rm in ( b ) 




FIG. 5. R/Rmin as functions of r in the modified Type-I f{R) models with n = 2, where (a), (b), 
(c), and (d) represent {gi,T) = (1,50), (10-2,50), (1,100), and (lO'^lOO), respectively. 

In Fig. O we display R/Rmin as functions of r in the modified Type-I f{R) models with 
n = 2, where Rmin = i^'^Tq{1 + t/tch)- The large values of gi and Tch correspond to the 
smaller amplitude of R/Rmin, and vise versa. If gi is small enough (see Figs. [5]d and[5li), 
the divergent behavior of the original viable f{R) part is more efficiency, and the amplitude 
increases until the effect of the R^ part is big enough. In Fig. El we illustrate the similar 
results as Fig.[5]for the modified Type-II f{R) models with two different values of gu, related 
to the background density. From the figures, we find that the amplitude in Fig. increases 
as that in Fig. [5^, and it is almost a constant in Fig. [6^. It is clear that the amplitude in 
the Type-II models is more sensitive than the Type-I ones when we change gj. 

In the physical system, we can estimate the order of gjj in different background densities in 
the modified Type-II f{R) models. Under the natural selection of the mass factor M ~ Mpi, 
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1.6| r- 



(b) 



1.5- 
1.4- 
1.3- 
1.2- 



0.9- 
0.8- 
0.7- 
0.6- 



\J 




T xlO"* T xlO"* 

FIG. 6. R/Rmin as functions of r in the modified Type-II f{R) models for Tch = 10"^ and {gjj, 13)- 
(1,25) and (10~2,20.4) for (a) and (b), respectively. 

leading to 



P = K'To/Rch > In 



Rii 



280 



(3.5) 



we find that g = 0(1), where we have assumed A ~ 0(1) and Rch/Mpi ~ A/M^ 
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This is a very low energy density since the physical density, such as that in the inner galaxy 
regime, is much larger. As a result, it is hard to detect the oscillation behavior and scalar 
mode of gravitational waves in the Type-II f{R) models. On the other hand, the curvature 
oscillation and scalar mode of gravitational waves could still be detected in the modified 
Type-I f{R) models. Note that the amplitude depends on both gj and r^h (see Figs. [5^ and 
[Sb). The galaxy collision takes billions years {tch ~ lO^^sec and Tch ~ 10^"^ ^ 1 for n = 2 
and P = 10^), so that 6R/Rmin when gj > 0(1). The curvature oscillation and scalar 
mode of gravitational waves can only be observed for a small value of n as it corresponds 
to gi < 0(1), which can be estimated by using Eq. (13. 3p with R^h ~ A ~ 10^'^^gcm^^. For 
example, for the inner galaxy (sun) with the density pm — 10'"^^ g/cm^{pQ ~ lAg/cm^), one 
finds that n < 11(= 1). 

In the above discussions, we have assumed the fiat space-time. Now, we would check 
whether this approach is valid or not when the curvature becomes large. If the high density 
region is described by the locally homogeneous and isotropic FRW metric, one has the 
form of g^i, = diag(— 1, a^, a^, a^), where a is the scale factor. Then, the gravitational field 
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(Friedmann) equation is given by 

3 (1 + fj,) = /tVm + I Wr - /) - 3HfR, (3.6) 
-2 (1 + fR)H = {pm - Pm) + fR- Hfn, (3.7) 

with R = 6 (^H + 2H^^, where H = d/a is the Hubble constant. In Eq. (13. 6p . the second 
term of is approximately equal to a constant curvature of XRj, which is much smaller than 
the first term (/t^p^), while the third term is 



As the large curvature oscillation corresponds to a small gi, we can estimate each quantity 
in Eq. (13.81) . Since Eq. (13.21) approximately describes a simple harmonic oscillator gjz" + 
z — {1 + T /Teh) = for (7/2^""^^ > 1, where z = R/k'^Tq, leading to the oscillating frequency 
00 = gj^^'^- Without loss of generality, we take g^z^^^'^ ~ 1. Consequently, the order of the 
oscillation amplitude of H'^ in Eq. ^B) is -?>HfR ~ -{H / Mpi)R{K^To)Ri <^ Rj <^ n^pm- 
Similarly, in Eq. (13. 7p . the oscillation amplitude of H is dominated by fji, which has the 
same order as the curvature oscillation amplitude: 

JL A = //^Y"^7^\ _R_ 

As a result, the singularity behavior comes from the small amplitude and high frequency 
of the scale factor oscillation. Hence, the result still holds under in the flat space-time 
limit. Explicitly, the covariant derivative yields the same result as the partial derivative, 
e.g., UR = -(r + 3HR^ ~ -R. 

Finally, it is worth to mention that the amplitude of R/Rmin would be converge if Tch < 
in the original Type-I f{R) models without the R^ term as shown in Fig. [71 This result 
can be seen as the time reversion of Tch > 0. This kind of singularity does not exist in the 
cosmological evolution. 

IV. CONCLUSIONS 

We have demonstrated that the singularity appears in the finite time in the viable f{R) 
gravity models. However, this singularity could be avoided by adding the /2"/M^'^"~^) term 
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FIG. 7. R/Rmin Vs. T when Tc/i = 
f{R) models without the R^ term. 



-1100 is a negative number with n = 2 in the original Type-I 



with 1 < n < 2 into the Lagrangian. We have exphcitly shown the cases with the most 
reasonable parameters of = 2 and M = Mpi. Even though we can prevent the divergence, 
the oscillation behavior still exists. If the oscillating amplitude is large enough, it would be 
detected by cosmological observations. The oscillating effect in the different scale depends 
on the types of f{R) models. In particular, in the Type-I models the effect can be induced 
in inner galaxy regime, but it is hard to distinguish the Type-II models from GR. 



As the simplest example, the scalar mode of gravitational waves is a typical phenomenon, 
which is different from GR, and it has the same origin as the curvature oscillation. As shown, 
the Type-I viable f{R) models can give good sources of the scalar gravitational waves when 
77, < 11. However, if n is too large (the case of n > 12 leads to (7 ~ 1 with the inner galaxy 
regime) or the background is too dense (denser than the solar density), the amplitude stays 
in a stable small fluctuation regime and it cannot be the source of the scalar mode. This 
behavior does not appear in the Type-II viable f{R) models, because the factor gjj is too 
large when k^Tq > 280Rch- This is still a very low density regime and hence, the energy of 



the scalar mode of gravitational waves 25| cannot be emitted. 
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